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O ■ ABSTRACT 

o ■ 

We show that the giant flares of soft gamma ray repeaters [E ~ 10^"^ erg) can 
push the inner regions of a fall-back disk out to larger radii by radiation pres- 
OO I sure, while matter remains bound to the system for plausible parameters. The 

subsequent relaxation of this pushed-back matter can account for the observed 
. _ enhanced X-ray emission after the August 27^^ giant flare of SGR 1900+14. Based 

I on the results of our models, we estimate that the ratio of the fluences of the 



enhanced X-ray emissions to that of the preceding bursts remains constant for a 



I particular SGR with similar pre burst inner disk conditions, which is consistent 

^ ' with the four different burst observations of SGR 1900+14. 
O " 

^ • Subject headings: pulsars: individual (SGR 1900+14) — stars: neutron — X- 

I ' rays: bursts — accretion, accretion disks 



1. Introduction 

Soft gamma ray repeaters (SGRs) are neutron stars that emit short (< 1 s) and luminous 
(^10^^ erg s~^) soft gamma ray bursts in their active phases. The burst repetition time scales 
extend from a second to years (see Hurley 2000 for a review). In their quiescent states, they 
emit persistent X-rays at luminosities similar to those of anomalous X-ray pulsars (AXPs) 
(L ~ 10^^ — 10^^ erg s~^). The spin periods of both SGRs and AXPs are in a remarkably 
narrow range (P ~ 5 - 12 s) (see Mereghetti 2000 for a review of AXPs). Four SGRs (and 
one candidate) and six AXPs are known up to date. Some of them were reported to be 
associated with supernova remnants indicating that they are young objects. Recently, some 
AXPs also showed bursts similar to those of SGRs, which probably imply that they belong 
to the same class of objects. 
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Over the burst history of SGRs, two giant flares were exhibited by SGR 0526-66 (Mazets 
et al. 1979) and SGR 1900+14 (Hurley et al. 1999). These giant flares are characterized 
by an initial hard spike with a peak luminosity ~ 10^^ — 10^^ erg s~^ which lasts a fraction 
of a second. The hard spike is followed by an oscillating tail that decays in a few minutes. 
Assuming isotropic emission the fluence of the entire giant flare is about ~ 10*^^ ergs (Hurley 
et al. 1999; Feroci et al. 2001; Mazets et al. 1999). 

Magnetar models can explain the super-Eddington luminosities of the normal and the 
giant bursts of SGRs by the sudden release of the very high magnetic energies from inside 
the neutron stars (Thompson & Duncan 1995). In an alternative class of models, fall-back 
disks around young neutron stars can account for the period evolution of these systems, and 
in particular for the period clustering of SGRs and AXPs (Chatterjee, Hernquist, & Narayan 
2000; Alpar 2001). Thompson et al. (2000) argued that the high luminosity of a giant flare 
would excavate any accretion disk to a large radius (due to the radiation momentum) and 
rebuilding of the entire disk takes months to years; so that the enhancement and the decay 
of the persistent X-ray flux after the giant flare could not be related to any disk accretion 
phenomenon. 

The persistent X-ray emission from SGR 1900+14 was reported to increase by a factor 
~ 700 about 1000 s after the giant flare. The subsequent decay is a power law with an index 
~ 0.7 (Woods et al. 2001). This increase and decay in the persistent X-ray emission of the 
SGR 1900+14 is our main interest here. It was proposed that the enhanced X-ray emission 
is due to the cooling of the neutron star crust after being heated by the energy of the giant 
flare (Lyubarsky, Eichler, & Thompson 2002). Here we show by means of a numerical disk 
model that (i) the X-ray enhancement can be explained in terms of the viscous relaxation 
of a disk pushed back by the giant flare, (ii) the amount of disk matter pushed out while 
remaining bound corresponds to a plausible fraction of the flare energy. The origin of the 
giant flare, which is probably the release of the high magnetic energy inside the NS by an 
instability, is not addressed in our model. 

In the next section, we summarize the X-ray observations of SGR 1900+14 reveahng 
the large flux changes in the persistent X-ray emission following the August 27 giant flare. 
In Sec. 3, we present the details of the numerical disk models. The results of the model flts 
are discussed in Sec. 4. The conclusions are summarized in Sec. 5. 
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2. The X-Ray Data 

The X-ray (2-10 keV) flux data following the August 27 giant flare of SGR 1900+14 
shown in Figures 1-3 was taken from Woods et al. (2001). The flrst and the second data 
points (filled squares) are from RXTE/ASM measurements, and correspond to ~ 24 minutes 
and ~ 2 hours after the giant flare (Remillard et al. 1998; see Woods et al. 2001 for the 
measurements and the associated uncertainties). The two data points about 20 days after 
the giant flare are from the net source intensity measurements of BeppoSAX (flUed triangle) 
and ASCA (flUed circle) satellites. The remaining data points (crosses) are estimated from 
the pulsed intensity measurements by RXTE/PCA (Woods et al. 2001) as follows. Four 
BeppoSAX NFI observations of SGR 1900+14 (2000 March/Aprfl, 1998 September, and 
1997 May) give similar pulsed fractions (~ 0.1) despite the varying intensity, pulse proflle, 
and burst activity. In the light of these observations. Woods et al.(2001) estimated the to- 
tal source intensity from the pulsed intensity measurements by assuming a constant pulsed 
fraction (^rms ~ 0.11). There is a good agreement between these estimates, and the Bep- 
poSAX and ASCA source intensity measurements about 20 days after the giant flare (see 
Woods et al. 2001 for more details). The reported relative pulsed fraction changes along the 
X-ray tails following the other observed bursts of SGR 1900+14 do not exceed a factor ~ 2 
in the extreme case (Lenters et al. 2003). This does not affect the quality of our fits, but 
might require a small modification of the model parameters presented here. Keeping these 
uncertainties in mind we adopt for our numerical model the X-ray data set shown in Figures 
1-3, obtained by scaling with the pulsed fraction T = 0.1 where only the pulsed signal is 
observed. 



3. The Numerical Model 

Assuming isotropic emission, the total emitted energy during the giant fiare is ~ 10^^ 
ergs (Mazcts et al. 1999). A fraction of this emission is expected to be absorbed by the 
disk depending on the solid angle provided by the disk for the isotropic emission. For such 
a point-like emission at the center of the disk, the radiation pressure is expected to affect 
mostly the inner regions of the disk by pushing the inner disk matter to larger radii depending 
on the energy imparted to the disk matter. This leads to large density gradients at the inner 
rim of the disk immediately after the giant ffare. We test whether the consequent viscous 
evolution of the disk can reproduce the X-ray ffux data, consistently with the reported energy 
arguments of the giant ffare. 

In our model, we represent pushed-back inner disk matter, which we assumed to be 
formed by the radiation pressure of the giant ffare, by a Gaussian surface density distribution 
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T,{R, t — 0) — Smax exp — ( \^° )^ , representing the pile up, added to the inner edge, at 
i?o, of the extended disk profile for which we chose the form E = T,o{Ro/R). Eq is a constant 
much less than Emax) R is the radial distance from the center of the disk, and Rq is the 
initial radial position of the center of the Gaussian. This form of the extended disk is close 
to the surface density profile of a standard thin disk (Shakura & Sunyaev 1973). In addition 
to the post-flare radius Rq, Sq, the Gaussian width and the maximum initial surface density 
Sniax (at the center of the Gaussian) are the free parameters of our model. The disk's 
inner radius (where the subsequent inflow of the pushed-back matter will be stopped 
by the magnetic pressure) , and the outer disk radius -Rout ^^P^ constant throughout the 
calculations. A constant outer disk radius was chosen due to numerical reasons. Outer disk 
properties can only affect the inflow rate through the inner disk after several weeks or more 
in the absence of large surface density gradients at the outer disk regions. We use the one 
dimensional disk code described in Ertan & Alpar (2002), originally constructed to simulate 
the black hole soft X-ray transient accretion disks in outburst. 

For a Keplerian thin disk the mass and angular momentum conservation equations give 
a non-linear diffusion equation for the surface density 



R'/'^ii^^R'/'] 



1) 



dE _ 3 d 
'dt ~ RdR 

(Frank et al. 1992), where u is the kinematic viscosity which, together with the surface 
densities, can be related to the disk midplane temperatures Tq through 

-T^ - (2) 
3r ^ 8 ^ ' 

T = kS is the vertically integrated optical depth, and a is the Stefan-Boltzmann constant. 

For the viscosity we use the standard a prescription u = acsh (Shakura & Sunyaev 1973) 

where cs = kTc/ fim-p is the local sound speed, n the mean molecular weight, h = cs/Oj^ the 

pressure scale height of the disk, and fij^ the local Keplerian angular velocity of the disk. 

We use electron scattering opacities (^es — 0.4 cm^ g"^)- We chose ji — 0.6 and a = 0.1 

which is typical of the hot state viscosities in the disk models of dwarf novae and soft X-ray 

transients. 

By setting x — 2it!^/^ and S — xE, Eq.(l) can be written in a simple form 

dS 12 

We divide the disk into 400 equally spaced grid points in x. This provides a better spatial 
resolution for the inner disk in comparison to a model with the same number of grid points 
equally spaced in R. 
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For a thin disk, the total disk luminosity is L^^^]^ — GMM-^^/2R-^^, and most of this 
emission comes from the inner disk, characterized by a disk black-body spectrum. Here, 
Mjj^ is the mass inflow rate arriving at the disk inner radius -Rj^, and M is the mass of 
the neutron star (NS). We take M = IAMq throughout the calculations. The accretion 
luminosity from the NS surface, = GMM^/R^, determines the observed luminosity in the 
X-ray band. The evolution of M^^{t) in the disk will be reflected in the accretion luminosity 
from the NS surface, depending on the fraction of matter accreted, / = M^/Mjj^ where 
is the mass accretion rate onto the star. We present three model calculations corresponding 
to different / values (0.1, 0.5, 0.9). 

While the observed luminosity is expected to be powered by accretion onto the NS 
surface, the spectra during the enhanced X-ray emission of SGR 1900-1-14 can be fitted by 
a single power- law (Woods et al. 2001). A scattering source, e.g. a hot corona, around the 
inner disk can significantly change the spectrum emitted from the neutron star surface and 
from the disk black-body spectrum into a power-law spectrum by means of inverse Compton 
scatterings. If the source of the corona is fed by the thermal instabilities at the surface (or 
inner rim) of the disk then the total luminosity remains constant for a given matter inflow 
rate and inner disk radius, while the spectrum may be modified from the input spectrum. 
Comparison of spectral models for emission from the NS surface or the disk with the observed 
2 — 10 keV band data may be misleading. We take the observed luminosity to represent the 
total luminosity assuming that most of the X-ray flux from the source is emitted in the 
observation band (2-10 keV). For the model fits, we relate the model luminosities to the 
fiuxes by F^^^]^ ~ (-^disk i)/(47rd^) and ~ L*/(47rd^) where d — 14.5 kpc is the 
distance of the source (Vrba et al. 2000). We set cos i — 0.8 and neglected the small time 
delay for the matter to travel from to R*. 



4. Results and Discussion 

The disk parameters for the model curves presented in Figures 1-3 are given in Table 1. 
The lower and the upper model curves in the Figures correspond to the fiuxes originating from 
the inner disk and from the NS surface respectively with = 2(M^/R^,){R^^/M^j^)L^j^g-^ — 
2/(i?jj^/i?*)Ljjgj^. For each of the three different / values (0.1, 0.5, 0.9) » L^^^y^. 
Our models produce good fits to the wide range of /. For each mass accretion ratio /, the 
quiescent luminosity gives the mass infiow rate in the disk. The R:^^ values given in Table 1 
are estimated Alfven radii for these mass infiow rates, taking the dipole magnetic moment 
II — 10^° G cm^. These results strongly suggest a viscously evolving disk origin for the 
observed post burst X-ray enhancement, but do not constrain /. In this range of /, rough 
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estimates with the thin disk model give m/ ~ 19 — 20.5 at the peak of the hght curve, and 
mi ~ 26, similar to the upper limits in the quiescent phase (Vrba et al, 2000). The upper 
limits placed by the IR observations about 8 days after the giant flare, when the X-ray flux 
has decreased to about one percent of its peak level, are mj> 22.8 and itlk^^ 20.8 (Kaplan 
et al. 2002). The IR expected light curve during the X-ray enhancement and in quiescence 
will be presented in a separate work. 

The energy given to the disk by the giant flare could be written as 5E = l3EAt ~ (310"^^ 
ergs where f3 = f]-^ + j3e is the fraction of the total flare energy absorbed by the disk. Part 
of the inner disk matter heated by the energy j3e5E can escape from the system, while the 
remaining part is pushed back by I3-^5E staying bound and piling up at the inner rim of 
the disk. /3 is expected to be around ~ 27r(2ifjj^)i?jj^/47ri??^ = H-^^/ R-^^ ~ few xlO~^ for a 
thin disk with M ~ lo^^-ie g g-i ^^ere i/j^ is the semi-thickness of the disk at R-^^. This 
ratio is roughly constant throughout the disk (e.g. Prank et al. 1992). The energy imparted 
by the flare to push back the inner disk matter is: SEy^ ~ {GM5M/2R^^\i — {R^^/Rq)\. 
This is almost equal to the binding energy, since we find that R^^/Rq ~ 1/3 for the models 
given in Table 1. The energy used up pushing back the disk is a fraction of the estimated 
energy, absorbed by the disk, f]-^-^ < j3. It is in fact likely that a larger amount of matter 
escapes from the system, than the amount 5M that is pushed back but remains bound, with 
/?e ~ (5 - 25)/5b. 

The maximum amount of mass that can escape from the inner disk during a burst can 
be estimated as ^Mj^gg ~ (2R-^.^/ GM)I35E ~ 10^^ g R-^^ g {13 /IQ^^) where R-^^ g is the inner 
disk radius in units of 10^ cm. During the lifetime of an SGR (~ 10^ yrs) which has a giant 
burst per century, the total mass loss would be 10^^ g R-^^ g (/3/10^^). 

If the pulsed fraction remains the same (~ 0.1) throughout the enhanced X-ray flux 
phase as estimated by Woods at al. (2001) we expect a connection between the mass inflow 
rate and the pulsed X-ray emission. In our models, the luminosity from the NS surface 
dominates the disk luminosity, and the pulsed fraction ~ 0.1 could be explained as the 
ratio of the emission beamed by the mass flow geometry through the polar caps to the 
isotropic emission from or near the NS surface. 

The time evolution in our models is quite prompt, with a viscous time scale ty ~ R? /v ~ 
10^ s, in agreement with the observed X-ray enhancement. Thompson et al. (2000) estimate 
a viscous time scale of ~ 10 yrs for the reestablishment of the inner disk mainly because they 
use pre burst mass flow rate M ~ 10^^ g s~^ in their estimate, instead of the appropriate post 
burst M, which is three orders of magnitude higher. Thompson et al. also take a — 0.01 and 
estimate the post burst inner disk radius to be Rq = 10^° cm. In our calculations, a = 0.1, 
typical of the outburst (hot) states of the soft X-ray transient and dwarf nova disk models. 
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The post burst pile up position i?o ~ 10^ cm in our models corresponds to the short viscous 
time scale. For smaller burst energies (IC^^^^^ ergs), the inner disk matter is pushed out to 
correspondingly smaller radii Rq, and t,^ could be as small as a few seconds. 

The enhanced mass inflow rate can modify the spin evolution significantly especially 
around the peak of the X-ray light curve when Lx is about ~ 700 times higher than its 
quiescent level. A similar increase in Mj^ corresponds to a decrease in the Alfven radius 
by a factor ~ 6 (or less if the field lines are compressed by rapid accretion). If the fastness 
parameter uj = Q^/Qj^{rji^) drops below unity a spin up is expected to prevail until uj has 
reduced back. If the system does not enter the spin-up phase or remains mostly in the spin- 
down phase in the high mass inflow regime, a sharp increase in the spin down torque could 
abruptly reduce the spin frequency. The changes in the spin evolution depends not only on 
the variations of the Mj^ and cu but also on how close the system is to rotational equilibrium 
in quiescence and on the inner disk structure (likely to be somewhat different from a thin 
disk geometry) during the X-ray enhancement phase. Due to these uncertainties it is hard 
to make a reliable estimate of the spin evolution during the short term unsteady phase of 
these systems. A detailed examination of the possible post-burst spin evolution considering 
different accretion geometries will be addressed in future work. 

Observations of four bursts from SGR 1900-1-14, including the August 27 giant flare 
and 3 smaller events, extending three orders of magnitude in flare fluence were studied by 
Lcntcrs ct al. (2003, sec especially their Fig. 13). Their study reveals that the ratio of the 
fluence of the enhanced X-ray emission 6Ex to the fluence of the preceding burst energy 
'^-^burst ^ 0.02, and remains constant from burst to burst. In our models, this ratio can 
be written as 7 = SEx/SE^^^.^^ ~ 2(3-^f{R^j^/R^) where both (3-^ and represent the pre 
burst inner disk conditions. (3-^ depends on the disk geometry and is very likely to be similar 
prior to the different bursts of SGR 1900-1-14. Our models with a constant / along the X-ray 
enhancement phase fits well to the data indicating that / remains constant along this phase. 
Since the X-ray enhancements following the three other smaller events trace accretion rates 
that were encountered along the decaying tail of the post giant fiare enhancement, a similar 
/ must be operating throughout the smaller enhancements following the three events. The 
remaining variable depends on the pre burst Mj^. An order of magnitude change in Mjj^ 
causes a change in R^^ by a factor < 2. So, based on our model results, it is understandable 
that the ratio 7 remains constant within a factor ~ 2 for different bursts of a particular 
SGR, consistent with the observations of SGR 1900-1-14. 7 may vary from source to source 
depending on the preburst inner disk conditions. 
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5. Conclusion 

We have shown that the X-ray flux curve following the 1998 August 27 giant flare of 
SGR 1900+14 can be accounted for by the enhanced accretion onto the neutron star surface 
due to the relaxation of the disk, starting from new initial conditions with the inner disk 
pushed back by a plausible fraction of the flare energy. For our disk models, the ratio of the 
fluence of the X-ray enhancement to the preceding burst energy remains roughly constant 
for bursts of a given SGR with similar preburst mass inflow rates, in agreement with the 
burst and enhancement observations of SGR 1900+14 (Lenters et al. 2003). This ratio can 
vary for different SGRs indicating their different inner disk conditions. 
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comments. We acknowledge support from the High Energy Astrophysics Research Group 
TBAG-(^-4 of TUBITAK and from the Astrophysics and Space Forum at Sabanci University. 
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Table 1: Model parameters for the flux evolution presented in Fig. 1. In all model calcu- 
lations, the viscosity parameter a = 0.1, the mean molecular weight fi = 0.6 and electron 
scattering opacities are used 





MODEL 1 


MODEL 2 


MODEL 3 


Smax (g cm 2) 


9.6 X 10^ 


3.0 X 10^ 


2.1 X 10^ 


Gaussian width (cm) 


2.4 X 10^ 


2.2 X 10^ 


2.2 X 10^ 


So/Smax 


0.012 


0.020 


0.022 


-Ro (cm) 


1.8 X 10*^ 


1.1 X 10*^ 


9.4 X 10* 


(cm) 


6.0 X 10* 


4.0 X 10* 


3.0 X 10* 


^out (cm) 


1.0 X 10^^ 


1.0 X 10^^ 


1.0 X 10^1 


/ 


0.1 


0.5 


0.9 


estimated j3-^ 


2 X 10'* 


5 X 10-5 


4 X 10-5 


SM (g) 


2 X 10^3 


3.5 X 10^2 


2 X 10^2 



